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ABSTRACT 

We present iV-body simulations of young substructured star clusters undergoing vari- 
ous dynamical evolutionary scenarios and examine the direct effects of interactions in 
the cluster on planetary systems. We model clusters initially in cool collapse, in virial 
equilibrium and expanding, and place a 1 Jupiter-mass planet at either 5 au or 30 au 
from their host stars, with zero eccentricity. We find that after lOMyr '^10 per cent 
of planets initially orbiting at 30 au have been liberated from their parent star and 
form a population of free-floating planets. A small number of these planets are cap- 
tured by other stars. A further ~10 per cent have their orbital eccentricity (and less 
often their semi-major axis) significantly altered. For planets originally at 5 au the 
fractions are a factor of 2 lower. The change in eccentricity is often accompanied by 
a change in orbital inclination which may lead to additional dynamical perturbations 
in planetary systems with multiple planets. The fraction of liberated and disrupted 
planetary systems is highest for subvirial clusters, but virial and supervirial clusters 
also dynamically process planetary systems, due to interactions in the substructure. 

Of the planets that become free-floating, those that remain observationally as- 
sociated with the cluster (i.e. within two half-mass radii of the cluster centre) have 
a similar velocity distribution to the entire star cluster, irrespective of whether they 
were on a 5 au or a 30 au orbit, with median velocities typically '^1 kms""'^. Conversely, 
those planets that are no longer associated with the cluster have similar velocities to 
the non-associated stars if they were originally at 5au ('^9kms~^), whereas the plan- 
ets originally at 30 au have much lower velocities (3.8kms~^) than the non-associated 
stars (lO.Skms^^). These findings highlight potential pitfalls of concluding that (a) 
planets with similar velocities to the cluster stars represent the very low-mass end of 
the IMF, and (b) planets on the periphery of a cluster with very difi'erent observed 
velocities form through difi'erent mechanisms. 

Key words: stars: formation - kinematics and dynamics - open clusters and associ- 
ations: general - planet-star interactions - planetary systems - methods: numerical 



1 INTRODUCTION 

A large pr oportion of stars form in clu stered environ- 
ments (e.g. iLada fc Ladall2003l : lLadall2010l ). Whether such 
star-forming regions are dense enough to undergo sig- 
nificant dynamical processing is currently under debate 
iBressert et al.ll2010l ). but dynamical considerations do sug- 
gest that some cluster s experience a dense phase dur- 
ing their evolution (e.g. iKr oupa 1995; Kroupa et al." 1999*; 
iMoraux et al. 200'7l: lAllison e t al, 2009; Parker et al. ,2009. ). 
Recentlv. lAllison et al.l ( 20091 ) showed that the mass segrega- 
tion in the Orion Nebula Cluster (ONC) can be explained if 
the cluster was originally subvirial, and substructured. This 
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causes the cluster to collapse in the first 1 Myr which leads to 
dynamical mass segregation an d also heavily processes the 
prim ordial binary population (jParker. Goodwin fc AllisonI 

Iml). 

Such a dense phase would have serious implications 
for the survivability of planetary systems: planets could 
either be directly liberated from their host stars after a 
close encounter with another star, or the orbital elements of 
the planets could be significantly alter ed leaving the \Aa.ne- 
tary s ystem dynamically unstable (e.g. Holman et al. 1993; 



Innanen et al. 1997; Takeda fc Rasio 2005; Malmbcrg et al.l 



20071 : iFabrvckv fc Tremaindl2007l : lNaoz et al.ll201ll ). 



In terms of density, the most extreme galactic stel- 
lar environments are Globular clusters, and A'^-body cal- 
culations have demonstrated that planets orbiting at dis- 
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tances of between 0.5 and 50 au woul d be liberated 
from their host stars in these clusters (|Hurlev fc Sharal 
l2002h . On the other hand, open clusters are not as 
dense but, as discussed above, may undergo a dense 
phase in which for example, the primordial stellar binary 
popu l ation can be signifi cantly processed jKroupa et alJ 
1 19991 : [ Parker et alJ l201lh . Several authors have stud- 



ied the effects of open cluster envi ronments on plane- 
tary s ystems fe.g. iLaughlin &: Adams! [l998: Bonnell et alJ 



200ll: Smith & Bonnelll 200ll: Davies & SiEurdsson 


2001; 


Adams et al.ll2006l; Freeeau et al.ll2006l: Sourzem et alJ 


200C, 



and references therein). In the main, these authors deter- 
mined the cross-sectional probability for scattering planets, 
based on the likelihood of an interaction with a passing sin- 
gle or binary star. 



A small number of authors performe d direct A'^-body 
simu l ations of clusters con taining planets jHurlev fc Sharal 
|2002| : ISpurzem et al.l l2009f ) , but mainly considered Globu- 
lar clusters. Almost all these papers assume that the clus- 
ter is in virial equilibriu m and has a sm ooth radial profile 

i e.g. a Plum mer sphere fPlummer ' 1 9 1 1') , or King profile, 
Kinelll96d ')') at birth. Adams ct al. (2006) considered sub- 
virial clusters, but with smooth Plummer spheres, rather 
than substructured environments. Observations of young 
star forming regions show them to be highly subst ructured 
(e.g. ICartwright fc Whi tworth"2004': Schmcia 2011) and of- 
ten subvirial fe.g. lPeretto et al.„2006 : Proszkow et al. 2QQ%. 

In this paper we focus on the dynamical evolution of 
substructured open clusters, with subvirial, virial and su- 
pervirial initial conditions. We perform A^-body simulations 
of such clusters and include planets directly in the models. 
We outline our method for setting up the clusters in Sec- 
tion O we present our results in Section and we discuss 
them in Section (4] placing them in the context of searches 
for free-floating planets and planets in open clusters and the 
fleld. We draw our conclusions in Section [S] 



2 METHOD 

In this section we describe our method of setting up star 
clusters with substructure, and the assignment of compan- 
ions (either stellar or planetary) to each star. 

2.1 Cluster set-up 

Observations of young star forming regions indi- 
cate that a high leve l of s ubstructure is present 
jCartwright fc Whitworthl 12004 ISanchez fc Alfard l2009l : 
ISchmeial20l'lf r A convenient way of creating substructure on 
all sc ales is the fractal distribution (iGoodwin fc Whitworthl 
|2004| ). which makes each location in the cluster statistically 
identical to any other. Note that we are not claiming that 
the best approximation of substructure in star clusters 
is the fractal method; rather that the fractal is the most 
convenient method of producing substructure. Its main 
advantage is that the substructure is deflned by just one 
number; the fractal dimension D. This defines how 'fractal' 
the cluster is, with D = 1.6 describing a highly clumpy 
cluster (in three dimensions) and D — 3.0 describing a 
uniform sphere. 

We set up the fractals according to the method in 



iGoodwin fc Whitworthl l|2004 ). This be gins by defining a 
cube of side Adiv (we adopt A'div = 2.0 throughout), in- 
side of which the fractal is built. A first-generation parent 
is placed at the centre of the cube, which then spawns A|i^, 
sub-cubes, each containing a first generation child at its cen- 
tre. The fractal is then built by determining which of the 
children will themselves become parents, and spawn their 
own offspring. This is determined by the fractal dimension, 
D, where the probability that the child becomes a parent is 
given by N^^~^\ For a lower fractal dimension fewer chil- 
dren will mature and the final distribution will contain more 
substructure. The mean number of maturing children is 2^, 
and so the preferred values of D are 1.6, 2.0, 2.6 or 3.0, which 
would correspond to the number of maturing children being 
an integer, and therefore res ult in fewer departures from th e 
specified fractal dimension (jGoodwin fc Whitworthll2004l '). 
Any children that do not become parents in a given step are 
removed, along with all of their parents. A small amount of 
noise is then added to the positions of the remaining chil- 
dren, preventing the cluster from having a gridded appear- 
ance and the children become parents of the next genera- 
tion. Each new parent then spawns A'^iv second-generation 
children in A^;^ sub-subcubes, with each second-generation 
child having a N^f^~^^ probability of becoming a second gen- 
eration parent. This process is repeated until there are sub- 
stantially more children than are required. The children are 
pruned to produce a sphere from the cube and are then ran- 
domly removed (so maintaining the fractal dimension) until 
the required number of children is left. These children then 
become stars in the cluster. 

We set up clusters with just one fractal dimension, 
D = 2.0, which gives the cluster a moderate level of sub- 
structure, but is not as extreme as, for example, a cluster 
with D = 1.6. The effect of varying the fractal dimension is 
to change the level of dynamical interactions tha t take place 
as the cluster evolves. lAUison et all l|2009l . I2OI0I ) show that 
the lower the fractal dimension, the more likely (and quickly) 
it is that dynamical mass segregation will occur. Trapezium 
systems can also form dynamically within 1 Myr in a clus- 
ter under going cool collapse with a low fractal dimension 
{D s; 2.0. 1 Allison fc Goodwinll201ll '). A higher proportion of 
primordial binary s ystems are also di srupted if the fractal 
dimension is lower (jParker et al.ll201lh . However, we adopt 
a mid-range value oi D = 2.0 throughout this work. 

To determine the velocity structure of the cluster, chil- 
dren inherit their parent's velocity plus a random component 
that decreases with each generation of the fractal. The chil- 
dren of the first generation are given random velocity com- 
ponents from a Gaussian of mean zero. The random com- 
ponent added to the children's velocity is also drawn from 
a Gaussian, but is then multiplied by 1/Adiv for each gen- 
eration of the fractal. This results in a velocity structure in 
which nearby stars have similar velocities, but distant stars 
can have very different velocities. The velocity of every star 
is scaled to obtain the desired virial ratio of the cluster. 

We examine three different evolutionary scenarios for 
our clusters. In the main, we ado pt the cool-co llapse scenario 
as advocated for the ONC bv lAUison et al.l 1^09, 201(J). 
Such clusters have a virial ratio of Q = 0.3, where we define 
the virial ratio as Q = (T and are the total ki- 

netic energy and total potential energy of the stars/planets, 
respectively). Therefore, a cluster with Q — 0.5 is in virial 
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Figure 1. The initial velocity distribution for our 'default' clus- 
ter, a subvirial {Q = 0.3), A^obj = 1500 fractal cluster. We 
have used the binary centre-of-mass velocity in the histogram, 
rather than the individual velocities. We show the median veloc- 
ity (0.98kms~^) in the cluster by the dot— dashed line. 



equilibrium and a cluster with Q — 0.3 is said to be sub- 
virial, or 'cool'. We also set up supervirial ('warm') clusters 
with a virial ratio of Q = 0.7. 

The initial velocity distribution for a subvirial {Q = 
0.3), A''obj = 1500 cluster is shown in Fig. \T\ We bin the 
centre-of-mass velocities, rather than the individual veloc- 
ities of the binary (either star-star or star-planet) compo- 
nents, as these represent the system velocities. The median 
cluster velocity (0.98 kms~^) is shown by the dot-dashed 
line. 



2.2 System properties 

Our 'typical' star cluster has a total o f 750 primary stars, 
with masses drawn from a 3-part (e.g. lKroupall2002l ) IMF 
of the form: 



N{M) 



M 
M 
M 



-0.3 
-1.3 
-2.3 



mo < M /Mq ^ mi , 
mi < M/Mq m2 , 
m2 < M/Mq s; ms , 



(1) 



and we choose mo = 0.08 M©, mi = 0.1 Mq, m2 = 0.5 Mq, 
and ms — 50 Mq. We adopt an upper limit of 50 Mq be- 
cause our default model is an Orion-like cluster, where the 
most massive star, 6^ Or i C, has a system mass of 45-50 Mq 
llKraus et "ai] |2007i . l2009'fl. Also, we randomly sample from 
the IMF, which in principle could lead to the most massive 
star dominating the mass of the e ntire cluster if we adopt 
an upper Umit of eith er 150 Mq (|Figeij |2005| ). or 300 Mq 
ICrowther et al.ll2010l ). In some simulations, we vary the 
number of stars, with one suite of clusters containing 1500 
primary stars initially, and another containing 100 primary 
stars. This samples two different cluster mass regimes, and 



^ 9^ OriC is likely to be a triple system, wh ere the most mas- 
sive component has a mass of 30 - 35 Mq l lKraus et al.ll2007l : 
iLehmann et al.ll2010l) . 



combined with our default model of 750 stars gives us three 
different cluster densities (as all the clusters have a radius 
of ~ 1 pc). 

Clusters are observed to have a —2 power law mass 
distribution from 10^ to 10^ Mq, so an equal mass of stars 
make up the Galactic field from all clusters (|Lada fc Ladal 
|2003). Our models therefore cover the first three orders of 
magnitude in cluster mass; we do not simulate more massive 
clusters due to computational limitations. 



2.2.1 Stellar companions 

For each primary star, we assign a companion based on the 
binary fraction associated with its mass. We divide primaries 
into four groups, roughly corresponding to the binary frac- 
tion of these systems observed in the Galactic fielcjf|. Primary 
masses in the range 0.08 ^ M /Mq < 0.47 are M-dw arfs, 
with a binary fraction of 0.42 (|Fischer fc Marcvl HoO^ ) . K- 
dwarfs have masses in the ran ge 0.47 ^ M/M q < 0.84 
and binary fraction of 0.45 ( Mayor et al.|[l993 ). We com- 
bine G-, F- and A-stars together, and so primary stars with 
masses from 0.8 4 ^ M/Mq ^ 2.5 are as signed a binary 
fraction of 0.57 (|Duauennov fc Mavor|[l99ll ). All stars more 
massive than 2.5 Mq are grouped together and assigned 
a binary fraction of unity, as massive stars ha ve a much 
large r binary fraction th a n low-mass stars (e. g. Abt et al] 
19901: iMason et all 1 19981: iKouwenhoven et al.l |2005| . l2007l : 



Pfalzner fc Olczakli2007riMason et al.ll2009l . and references 



therein) . 

Secondary masses in these stellar binary systems are 
drawn from a flat mass ratio distribution; recent work by 
iReggiani fc Meveil (|201ll ) has shown the companion mass 
ratio of binary stars in the field to be consistent with being 
drawn from a flat distribution, rather than random pairing 
from the IMF. We note that drawing companions from a flat 
distribution means we do not recover a Kroupa IMF. 

In accordance wit h t he observatio ns of 

iDuauennov fc Mavoj l|l99ll ) and iRaghavan et all |20lO), 
the periods of binary star systems are drawn from the 
following period generating function: 



/ (logioP) cx exp 



-(logioP - logioP) 
2(t2 



(2) 



where logioP — 4.8, aiogjoP = 2.3 and P is in days. We 
convert the periods to semi-major axes using the the masses 
of the two components (with logioP = 4.8 corresponding to 
a semi-major axis of roughly 30 au). 

The eccentricities of intermediate- and wide-separation 
stellar binaries in the field are well approxim ated by a ther- 
mal distribution (|Heggidll975l : lKroupall20oi ): 

Me) = 2e. (3) 

In the sample of IDuauennov fc Mayor! l|l99ll ). short- 
separation binaries are observed to be on circular orbits, 
which we account for by reselecting the eccentricity of a 

^ We note that the binary fraction in most cluster s may be much 
higher, in som e cases approaching 100 per cent l lKroupal |2008| : 
lGoodwi"nll2010l . and references therein). We will address the issue 
of planets orbiting the component (s) of stellar binary systems in 
star clusters in a future paper. 
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system if in exceeds the following period-dependent value 
(Parker fc Goodwinllioogl ): 



etid = - [0.95 + tanh (0.61ogioP - 



2.2.2 Planetary companions 



1.7)] . 



(4) 



Stars which do not have a stellar companion based on the 
above criteria are then all assigned a planetary-mass com- 
panion. These planets are all 1 Jupiter mass (IMjup = 
9.5 X 1O~*M0), are placed on circular orbits (e = 0) and have 
the same initial semi-major axes. In certain simulations, the 
planets are placed at 5 au, corresponding to a Jupiter-like 
orbit, or 30 au, corresponding to a Neptune-like orbit. Due 
to the stellar system constraints above, no planets are placed 
around stars with masses exceeding 2.5 M©. 

Comparing our model setup to detected exoplanetary 
systems in terms of gravitational forces between the planet 
and the host star we note the following: a 1 Mjup object 
around a solar mass star at 5 au is similar to t he directly 
imaged exoplanet HR8799 e (|Marois et al.ll20l"ol ') assuming 
its current projected separation corresponds to its true semi- 
major axis. For a 30 au orbit the gravitational force is com- 
parab le to the 'b' planet of the HR8799 system (|Marois et al.l 
I2OO8I ). The planets c and d are intermediate case Other 
directly imaged planets are eit her much more wea kly bound, 
as is the case for Fomalhaut bjKalas et al.ll2008l ) and IRXS 



J 1609 b (jLafreniere et alj 20101). or more tightly bound, 
the C Eise of /? Pictoris b ( Lagrange et al.ll2010l : iQuanz et al.l 
|2010| ). compared to our model setup. 

Whilst creating planetary systems with only one planet 
is a very simplistic model, creating systems with more than 
one planet would prohibitively increase the run-time of our 
simulations, and we consider any change in orbital parame- 
ters of a Jupiter-mass planet to be indicative of the general 
effect of encounters in a star cluster on a fully populated 
planetary system. 

The simulations are run using the kira integrator in the 
Starlab package (e.g. iPortegies Zwart et af]|l999l . bOOll ) for 
10 Myr. We do not include the effects of stellar evolution in 
the simulations. A summary of all the simulations used in 
the paper is given in Table [1] 



3 RESULTS 

In this section we will focus on our 'default model'; a sub- 
virial (Q — 0.3), A'' = 1500 object cluster. We will describe 
the dynamical evolution of this cluster before describing the 
effects of this evolution on planetary systems. We will then 
outline the differences in the number of affected planetary 
systems for the different initial conditions. 



3.1 Cluster evolution 

In Fig. [2] we show the morphology of a typical cluster un- 
dergoing cool collapse. In panel (a) we show the cluster be- 



^ These estimates assume an age of 30 Myr for the HR8799 sys- 
tem in which case the masses are roughly 7 Mjup for the planets 



Table 1. A summary of the different cluster properties adopted 
for the simulations. The values in the columns are: the number 
of stars and planets in each cluster (A^obj): the typical mass of 
this cluster {Mt.iustcr)i the virial ratio (and corresponding virial 
'state'; 'cool' (c), virialised (v) or 'warm' (w)), and the separations 
of the planets. 





-^cluster 


Q 


Planet separations 


1500 


~ 6 X 10^ M0 


0.3 (c) 


5 au 


1500 


~ 6 X IO^Mq 


0.3 (c) 


30 au 


1500 


~ 6 X 10^ M0 


0.5 (v) 


30 au 


1500 


~ 6 X 10^ Mq 


0.7 (w) 


30 au 


200 


~ 10^ M0 


0.3 (c) 


30 au 


3000 


~ 10^ M0 


0.3 (c) 


30 au 



fore dynamical evolution, and the substructure is clearly ev- 
ident. After 1 Myr (panel (b)) the substructure has almost 
been erased, although the cluster has a similar spatial ex- 
tent. However, after 10 Myr (the time at which we analyse 
the fraction of affected planetary systems) the cluster has 
expanded following the dense phase of the cool collapse. 

We show the evolution over 10 Myr of the half-mass ra- 
dius, ri/2 of the cluster in Fig. 3(a) 
Pccnt in Fig. 



3(b) 



and the central density. 
We calculate the half-mass radius from 
the centre of mass of the cluster, and we define the central 
density as 

0.5A/clustor 



Pccnt — 



4 3 



(5) 



d, and e, and 5 Mjup for planet b. 



where Mduster is the total mass of the cluster. 

As the cluster collapses, the half-mass radius reaches a 
minimum of 0.4 pc, whereas the central density of the cluster 
rapidly increases in the first Myr, peaking at 1200 Mq pc~^ 
before the cluster expands and relaxes. The density after 
10 Myr is 24M0pc~^. Although the virial and supervirial 
clusters in our analysis do not reach such high densities, 
we still expect some dynamical interactions within the sub- 
structur e, as found for clu sters with a high stellar binary 
fraction (|Parker et aLlbOllI ). 

As the cluster evolves, a fraction of planets and stars 
become unbound from the cluster. We show the fraction of 
systems that remain bound in Fig. U Firstly, we determine 
whether a planet or star is bound to the cluster if it has a 
negative binding energy with respect to the centre of mass 
and velocity of the clusteiQ. The fraction of such bound plan- 
ets is shown by the solid line in Fig. 2] and the fraction of 
bound stars is shown by the dashed line. Secondly, we cal- 
culate the fractions of planets and stars that reside within 
two half-mass radii of the centre of the cluster; an observer 
may not have enough information to determine whether an 
object is energetically bound to the cluster and this second 
constraint provides a conservative lower limit to the num- 
ber of 'bound' systems. The fractions of planets and stars 
within 2 ri/2 are shown in Fig.|4]by the dot-dashed, and dot- 

^ If the star or planet is in a binary system then we use the 
centre of mass velocity of the binary, as a planet, or secondary 
component of a binary may have an orbital velocity in excess of 
the cluster escape velocity. 
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Figure 2. Typical morpiiologies for tlie 'default' model (Q = 0.3, A^obj = 1500) clusters at (a) OMyr, (b) IMyr and (c) lOMyr. The 
cluster is initially substructured, but this is erased in the first Myr through dynamical interactions in the substructure and the global 
collapse of the cluster. Note the difference in spatial extent between 1 and 10 Myr. 




Time (Myr) 
(a) Half-mass radius 

Figure 3. Evolution of (a) the half-mass radius, S'^d the (b) central density, pccnt, 

cluster over the duration of the simulation. The cluster reaches its densest phase (pcent = 



Time (Myr) 
(b) Density 
of the 'default' model {Q = 0.3, N^i^j 



1500) 



1200 Mqpc"'^, ri/2 = 0.41 pc) at IMyr. 



ted lines respectively. Using this second criteria, we see that 
at 10 Myr only 70 per cent of stars and planets remain within 
2ri/2 of the cluster centre. We will refer to systems within 
2ri/2 of the cluster centre as 'associated' with the cluster, 
and systems outside of this radius as 'non-associated'. 

3.2 Liberated and disrupted systems 

In our simulations, we determine whether a star or planet 
is in a bo und system using th e nearest neighbour method 
outlined in lParker et al.l l|2009r) . If two stars, or a star and a 
planet, are mutual nearest neighbours, and have a negative 
binding energy, we deem them to be a bonafide binary sys- 
tem. We also track planets in triple systems; if the two stars 
and the planet, or star and two planets are all mutual nearest 
and second nearest neighbours, then they are an observed 



triple. Note that such systems are often transient and gen- 
erally hierarchical, where the inner separation Oin is much 
smaller than the outer separation, aout (dout >> oiin). When 
this criterion is fulfilled, we determine the binding energy of 
the inner orbit, and then the outer orbit (using the centre 
of mass of the inner system). If the outer binding energy is 
positive then we do not include this system as a triple in our 
analysis. 

In the following analysis, we will refer to a planetary 
system being 'disrupted' if it suffers a significant change in 
either eccentricity or separation (we define what we consider 
to be reasonable thresholds below) but remains bound to its 
host star, or 'liberated' if it is removed from its host star 
through dynamical interactions. 

Planets that have been liberated can subsequently be 
captured by a star that is not their parent. Additonally, 
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Figure 4. The fraction of systems bound to tlie cluster as a 
function of time. We show the fraction of planets and stars that 
are energetically bound to the cluster's centre of mass (the solid 
and dashed lines, respectively); and the fraction of planets and 
stars that are located within two half-mass radii of the cluster 
centre and therefore would likely be observationally categorised as 
cluster members (the dot-dashed and dotted lines, respectively). 



the interaction of two systems can result in planets being 
exchanged between stars. As these events are both relatively 
rare, we refer to planets being 'captured' if they experience 
either event. We discuss such systems in Section [3.31 

In Fig. [S] we show the effects of dynamical interactions 
on planetary orbits in clusters undergoing cool collapse. In 
this plot we have combined the results of 10 different clus- 
ters, identical apart from the random number seeds used to 
initialise the simulations. Th e reasons for this are two-fold: 
(i) we improve the statistics (|Adams et al.l (|2006l ) find that 
averaging 10 simulations together produces the required sta- 
tistical significance), and (ii) the field is the sum of differ- 
ent star formation regions, so once these clusters disperse 
it becomes difficult to differentiate between different forma- 
tion regions based on the observations of the field alone. In 
Fig. [5] we plot eccentricity against separation for systems 
containing planets after 10 Myr, for planets initially at 5 au 
and planets initially at 30 au (Fig. 



(Fig. 5(a) 



5(b) I. The 



open circles denote primordial planetary systems, where the 
planet is still orbiting its parent star. The crosses indicate 
captured systems, where the planet is not orbiting its parent 
star. Finally, the asterisks indicate planets in triple systems; 
the plotted semi- major axis can either be ain or Oout, de- 
pending on whether the triple consists of a star orbiting a 
star-planet binary, or a planet orbiting a star-star binary, 
respectively. In the rare instance of a triple containing 2 or 
3 planets, we plot the inner-most semi-major axis, as this is 
less susceptible to dynamical break-up. 

The orbital eccentricities can be excited from to al- 
most 1 in some cases, and in fewer systems the planetary 
separations are either hardened or softened. Comparing the 
systems at 5 au to those at 30 au, we see that the effect is 
much more pronounced for the systems at 30 au. 

In order to quantify the effects of cluster evolution on 
the planets in our simulations, we define two criteria for a 
planetary system to be 'disrupted': either the eccentricity is 



raised from zero to 0.1, or the semi- major axis of the planet 
decreases or increases by 10 per cent or more. The two pro- 
cesses are by no means mutually exclusive; a system with 
an altered semi-major axis is likely to have a non-zero ec- 
centricity. This eccentricity t hresh old was chosen in particu- 
lar because [Raymond et al.l (|201ll ) recently showed that the 
eccentricity of the innermost giant planet can have signifi- 
cant impact on the evolution of forming planetary systems 
and the chances of survival for terrestrial planets closer to 
the star. Specifically, these authors modeled the evolution of 
planetary systems in which the giant planets in the system 
are subject to dynamical instabilities, and they found that in 
the case of a final eccentricity ^0.1 for the innermost giant 
planet, the vast majority of the simulations ended up with 
unstable systems or systems containing only one terrestrial 
planet. It should be noted that this change in eccentricity 
would not cause a giant planet at 30 au to reach the same pe- 
riapsis as the terrestrial regime of its planetary system; the 
planet would require its eccentricity to be raised to e > 0.85 
for this to occur. 

Our threshold for the change in semi-major axis is 
rather arbitrary; however, as we will see, systems with their 
eccentricities raised to above 0.1 almost always have their 
semi-major axis altered by 10 per cent or more and would 
therefore be 'disrupted' based on our first criterion alone. 

Adding together the simulations in which we place the 
planets at 5 au initially, we find that of the 4123 birth plan- 
etary systems, 198 (5 ± 2 per cent) planets are liberated 
from their host star, one of which is captured by another 
star; whilst 3925 (95 ± 2 per cent) survive the 10 Myr of 
cluster evolution. Of these 3925 'preserved systemsQ 104 
(2.6 ± 1.3 per cent) have eccentricities excited above 0.1, 
and 56 (1.4 ± 1.0 per cent) have their separations altered by 
more than 10 per cent (±0.5 au). The number of systems 
that have their eccentricity and semi-major axis altered be- 
yond the thresholds is 55 (1.4 ± 1.0 per cent). 

As is readily apparent from inspection of Fig. [5] plan- 
ets on Neptune-like orbits (at 30 au) are more susceptible 
to disruption. In this scenario, of the 4148 birth planetary 
systems, 503 (12 ± 3 per cent) planets are liberated, sev- 
enteen of which are captured by another star, whilst 3645 
(88 ± 3 per cent) survive the 10 Myr of cluster evolution. Of 
these 3645 preserved systems, 398 (11 ±4 per cent) have ec- 
centricities excited above 0.1, and 189 (5 ± 2 per cent) have 
their separations altered by more than 10 per cent (±3au). 
187 (5 ± 2 per cent) systems have both their eccentricity and 
semi-major axis altered, suggesting that any system that has 
its semi-major axis altered will also experience a change in 
eccentricity, whereas a system can suffer a change in eccen- 
tricity but its semi-major axis will remain unaffected. The 
fraction of planets originally at 30 au that would directly 
penetrate the terrestrial regime (e ^ 0.85) is 0.8 ± 0.6 per 
cent. 

The inclination angle of a planet's orbit can also be 
significantly altered during cluster evolution, as shown in 
Fig. |S] For the planets originally at 30 au, we plot inclina- 
tion against eccentricity, and we see that even planets that 
do not undergo strong changes in eccentricity can have a 



^ Note that we are not considering systems that formed via cap- 
ture during the evolution of the cluster. 
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(a) Planets at 5 au (b) Planets at 30 au 

Figure 5. Orbital eccentricity versus separation for IMjup planets initially at (a) 5au, and (b) 30 au, for subvirial N^hj = 1500 clusters. 
We have summed together the results of 10 different clusters with the same initial conditions after 10 Myr of evolution. Circles indicate 
'birth' planetary systems, crosses show captured planets and the asterisks show triple systems containing planets. 



100 r 




Figure 6. Inclination in degrees versus eccentricity for all pre- 
served planetary systems (the circles) originally on 30 au orbits 
(the planets all have a birth inclination of zero), in our default 
^obj = 1500, Q = 0.3 cluster. We also show the small number of 
captured systems by the (red) crosses. The values shown are at 
time 10 Myr. 



non-zero inclinatiorffj. Most interestingly, planets which ex- 
perience a significant change in eccentricity (i.e. e > 0.1) or 
are captured tend to have the highest inclinations. We will 
return to this in Section 



^ We also see the same behaviour for the planets originally on 
5 au orbits. 



3.2.1 Dependence on initial virial ratio 

The results described above are for planets in clusters which 
evolve via the cool collapse scenario. At present, it is unclear 
what proportion of clusters undergo cool collapse, and the 
respective proportions of clusters that are supervirial and 
expand from an early age, and those that remain in virial 
equilibrium. We also conduct simulations in which the over- 
all virial ratio of the cluster was 0.5 (virial equilibrium) and 
0.7 ('warm'). For simplicity we only consider the clusters 
with planets initially at 30 au, and we compare them to the 
clusters in cool collapse (virial ratio of 0.3). The results are 
shown in Fig. [T] In Fig 



7(a; 



ated systems at 10 Myr, in Fig. 7(b) 



we show the fraction of liber- 
we show the fraction 



of preserved systems that have their eccentricities raised to 



above 0.1 and in Fig. 7(c) we show the fraction of preserved 



systems that have their separations altered by more than 10 
per cent. 

As we would expect, the clusters undergoing cool col- 
lapse (and therefore subject to a very dense phase) are more 
likely to disrupt planetary systems than clusters that are in 
virial equilibrium or expanding. One might expect that an 
expanding cluster would not process the planetary popu- 
lation at all; however whilst the cluster is expanding the 
stars in the subclumps in the fractal interact and decay on 
a timescale much less than the time taken for the cluster 
to expand. Therefore these systems undergo dynamical in- 
teractions, albeit at a lower level than the clusters in virial 
equilibrium or in cool collapse, as demonstrated in Fig. [T] 



3.2.2 Dependence on density 

The fraction of planetary systems that are liberated, and 
have their eccentricities and semi-major axes altered varies 
as a function of the density of the cluster. We plot the 
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Figure 7. The fraction of (a) hberated systems, (b) preserved systems with eccentricity > 0.1, and (c) preserved systems with semi-major 
axis altered by more than 10 per cent, as a function of the initial viriai ratio of the clusters. We show clusters with initial viriai ratio 
Q = 0.3 (subvirial, or 'cool', which collapse), Q = 0.5 (viriai equilibrium) and Q = 0.7 (supervirial, or 'warm', which expand). The 
clusters have A^obj = 1500 and the planets were originally on 30 au orbits. The values shown are at time 10 Myr. 
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Figure 8. The fraction of (a) liberated systems, (b) preserved systems with eccentricity > 0.1, and (c) preserved systems with semi-major 
axis altered by more than 10 per cent, as a function of the maximum density reached by the subvirial (Q = 0.3) clusters. This corresponds 
to a time of 0.5 Myr for clusters containing 3000 objects, 1 Myr for 1500 objects, and 2.5 Myr for the clusters with 200 objects. The values 
shown are the fractions of liberated/disrupted planets at 10 Myr. 



results (at 10 Myr) for clusters with Q = 0.3 and plan- 
ets originally at 30 au in Fig. [S] In each panel we show 
the values for clusters with (from left to right) 3000, 1500 
and 200 objects. The critical value to be considered here 
is the most dense phase of the cluster's evolution, which 
occurs at ~ 0.5 Myr for clusters containing 3000 objects, 
^ 1 Myr for clusters containing 1500 objects, and at ~ 
2.5 Myr for clusters containing only 200 objects. The clus- 
ters containing 3000 objects obtain maximum densities of 
~ 3420 M© pc~^ (~ 7800 stars pc"^), whereas at the other 
end of the scale the clusters containing 200 objects reach 
densities of ~ 330 Mq pc~^ ('^ 680 stars pc"''). In the most 
dense clusters, the fraction of systems with altered eccen- 
tricities is 18 per cent, compared to only 4 per cent for the 
least dense clusters. The fraction of systems with altered 
semi-major axes follows a similar trend, with 6 per cent of 
systems having altered separations in the most dense clus- 
ters, compared to only 2 per cent for the least dense clusters. 



3.3 Captured systems 

As mentioned in Section [3.21 there are several systems where 
planets liberated from their host stars are captured by an- 
other star (the crosses in Fig. [5|. Whilst there is only one 
such system for the planets originally at 5 au, there are sev- 
enteen in the 30 au case. All of these systems are 'disrupted' 
in the sense that their eccentricity is > 0.1 and their semi- 
major axis also deviates by more than 10 per cent from the 
initial value. However, and particularly in the 30 au case, 
they do not populate any specific region in the plot in Fig- 
ure [S] making it observationally impossible to distinguish 
between captured or just disrupted primordial planets from 
the orbital elements alone. 

There is however, a high fraction of captured planets 
on retrograde orbits. In Fig. [5] we show (by the crosses) the 
orbital inclinations of planets originally on 30 au orbits that 
have been captured by another star. Of these seventeen plan- 
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ets, 8 (42 per cent) are on retrograde orbits (> 90°). This 
compares to the fraction of all planets that are on retrograde 
orbits of only 1 per cent. 



3.4 Velocities of liberated planets 

In Fig. [9] we show the distribution of velocities for planets 
that are single after lOMyr, i.e. those that have been lib- 
erated from their host stars and are now free-floating. We 
determine whether a free-floating planet is still associated 
with the cluster (within two half-mass radii of the cluster 
centre - see Fig. |4)), and then plot the velocity distribu- 
tion of these single planets in Figs. 9(a) and 9(b) Fig. 9(a) 
shows the distribution of planets originally on 5 au orbits, 



and Fig. 9(b) shows the distribution of planets originally on 
30 au orbits. 

We show the median of the free-floating planets' ve- 
locity distribution by the dashed line, and the distribu- 
tion of all of the cluster members by the dot-dashed lin4l|. 
The planets originally on 5 au orbits have a median veloc- 
ity of 0.94 km s^^, compared to the median cluster velocity 
of 0.87 km s~^. Planets liberated from 30 au orbits have a 
median velocity of 0.78 kms~^, as does the entire cluster. 

We also compare the velocity distribution of free- 
floating planets that are no longer associated with the clus- 
ter (at a distance of more than two half-mass radii from the 
centre) to the velocity distribution of every object which 
is also no longer associated with the cluster. In Fig. 9(c) 
we show this distribution for planets originally on 5 au or- 
bits, and the corresponding plot for planets originally on 
30 au orbits in Fig. 9(d) We see that the median non- 
associated planet velocity (9.0 kms~^; the dashed line in the 
plots) is similar to the median non-associated object veloc- 
ity (8.5kms~^; the dashed line) for planets originally on 
5au orbits. The median velocity for non-associated planets 
originally on 30 au orbits is 3.8kms~^, whereas the median 
velocity for all non-associated objects is 10.8 km s~^. 

We will discuss these results further in the following 
Section. 



4 DISCUSSION 

From inspection, o ur results are quali t atively similar to 
those obtained by iLaughlin fc Adamd (| 19981 ). who per- 
formed cross-sectional scattering experiments on planetary 
systems assuming an orbital cross section for disruption of 
(cr) = (230au)^, a typical cluster density n of 1000 stars pc""^ 
and velocity dispersion n of 1 kms~^. They performed Monte 
Carlo scattering experiments on Jupiter-mass planets at 5 au 
with initially zero eccentricity, and noted a spread in or- 



bital parameters similar to our Fig. 5(a) However, our clus- 



ters can reach higher densities (up to 10* sta rs pc ^) when 



in cool collapse (Fig. [8 l ). Later authors (e.g . Bonnell et al.l 
I2OOII : ISmith fc Bonnelfl200ll : lAdams et aLri2006l ) assumed 



densities consistent with the Orion Nebular cluster (ONC) 
of 10'' — 10'* stars pc"'' but with smooth and mainly virialised 



^ For stars or planets in binary systems, we use the centre-of-mass 
velocity, as discussed in Section 2.1. 



initial conditions. We expand upon this earlier work by look- 
ing at the effects of dynamical evolution in substructured, 
and sub-virial, super-virial and virialised clusters. 

The majority of planets remain unscathed during the 
first 10 Myr of the cluster's evolution, with typically 90 per 
cent of planets at 30 au surviving break-up, and 95 per 
cent of planets at 5 au surviving. However, of these 'pre- 
served' planetary systems, a further 10 per cent would ex- 
pect to have their eccentricities raised by more than 0.1. 
The fraction of systems that have their separations altered 
by more than 10 per cent is typically a factor of two lower 
than the systems that have their eccentricities altered (see 
Figs. [7] and [HI). 

Whilst very few planets originally on 5au orbits are 
disrupted, there are a number that have high eccentrici- 
ties. Quite possibly, these planets will eventually be ejected 
from the system, as they are likely to cross the orbits 
of other planets in the system, leading to planet-planet 
scattering events (e.g. iMalmberg. Davies fc Chambersll2007l . 
and references therein). It is unlikely that inner terres- 
trial planets would be preserved or even able to form if a 
Jupiter-mass planet a t 5 au had a large orbital eccentricity 
l|Ravmond et al.ll2011^ . 

Of the larger fraction of planets originally at 30 au that 
have their separations and eccentricities increased, there 
is a significant population that occupy a phase space in 
which they might be directly detected in future imaging 
surveys carried out with dedicated planet-finding instru- 
ments such as Spectro-Polarimetric High-con trast Exoplanet 
REsearch (SPHERE) at the VLT (|Beuzit e t al. 2006) and 
Gemini Planet Imager (GPI) at the Gemini Observatory 
(jMacintosh et al1l2006l ). In particular the planets with high 
eccentricities and large semi-major axes spend a signifi- 
cant proportion of their orbital period far away from their 
host star. Surveys carried out in recent years are yet to 
reach the required sensitivity to det ect objects with masses 
below a few Jupite r masses (e.g., iLafreniere et al. 
IChauvin et al.l [20101 ). but this may change with the next 
generation of instruments. If a population of distant planets 
is detected and if formation mechanisms preclude the forma- 
tion in-situ of planets at such separations, then their origin 
could be dynamical. We note that a significant number of 
systems (25 per cent) with e > 0.3 and a > 50 au are actu- 
ally captured systems, implying that planets observed with 
these orbital characteristics may not be orbiting their parent 
star. This possibility should also be kept in mi nd when the 
origin of the distant exop lanets Fomalhaut b |Kalas et al.l 
1 200k ) and IRXS J 1609 b (|Lafi-eniere et al.ll2010l ') with pro- 
jected separations of ~115 au and ~330 au, respectively, is 
discussed. 

A large number of preserved systems have inclination 
angles raised due to interactions in the cluster (see Fig. [6]). 
Whilst the change in angle (from i = 0° at birth) for the ma- 
jority of these systems is only of order a few degrees, systems 
which we classify as being dynamically disrupted (based on 
e > 0.1) are likely to have inclination angles in the range 
where they could be lead to secondary dynamical processes, 
such as the Kozai mechanism (Kozai 1962t). The Kozai mech- 
anism can operate if the inclination angle of an outer body 
is raised to within the range 39.23° < ikoz < 140.77° of the 
orbital plane of an inner body. The total fraction of pre- 
served systems that have their angles excited to this range 
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Figure 9. Distribution of velocities of liberated, or free-floating, planets initially at (a) 5au, and (b) 30 au, which are still associated 
with the cluster — i.e. within two half-mass radii of the centre, for subvirial iVobj = 1500 clusters. We have summed together the results 
of 10 different clusters with the same initial conditions after 10 Myr of evolution. We also show the velocities of planets originally on 
(c) 5 au, and (d) 30 au orbits, that are no longer within two half-mass radii of the cluster centre. The median velocities of the liberated 
planets are shown by the dashed lines. The median velocities of all objects (either associated — panels (a) and (b); or non-associated - 
panels (c) and (d)) are shown by the dot— dashed lines. 



is only 2 per cent; however, for the systems with e > 0.1 the 
fraction is much higher, at 23 per cent. Additionally, of the 
captured systems, 76 per cent have inclination angles in the 

Kozai range. 

Recently, iNaoz et al.l l|201lh showed that Hot Jupiters, 
i.e., gas giant planets in a very close-in prograde or ret- 
rograde orbit, may result from secular planet-planet inter- 
actions. In their simulations a gas giant planet is initially 
orbiting on an almost circular orbit at several au, while a 
second gas giant planet is orbiting on a distant (a few tens 
of au), eccentric, and highly inclined orbit. Angular momen- 
tum exchange and eventually tidal dissipation may then lead 
to eccentricity fluctuations, orbital decay and circularisation 
of the first planet. In some cases even the orbital inclination 
was significantly altered leading to a retrograde orbital mo- 



tion. Our results show that dynamical interactions taking 
place in the cluster phase are one possibility to provide the 
initial conditions for this mechanism to worlfl. Therefore, in 
a more general sense, the altered eccentricities, semi-major 
axes and inclinations observed in our simulations may not be 
representative of the dynamical end-state of these planetary 
systems. 

The fraction of planetary systems that are affected by 
direct dynamical processing is dependent on density but also 



* An alternative mechanism to provide these initial conditions, 
where gas is captured onto a protoplanetary disc an d alters the 
inclin ation angle of the disc, was recently proposed bv lThies et al.l 
ll201ll) . 
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on th e init ial virial rati o , Q, of the system. lAllison et alj 
l|2009l ') and IParker et all (|201ll ') have shown that the 'cool 
collapse' scenario can explain the observed level of mass seg- 
regation, and the binary separation distribution in the ONC, 
if the cluster is initially substructured. If a significant pro- 
portion of clusters evolve in this fashion, then we would ex- 
pect a significant fraction (12 per cent) of free-floating plan- 
ets per cluster, and 10 per cent of the planets to have altered 
orbital parameters. This compares to around 10 per cent of 
extrasolar planets that could be affected by secondary dy- 
namical effects (such as the Kozai m echnism) in a dense 
cluster environment (jParker &: Goodw in 2009). If the clus- 
ter is initially in virial equilibrium, or is supervirial, fewer 
planetary systems are affected, although the clumps in the 
substructure dynamically decay in all cases, meaning that 
planets will be affected in almost all clustered environments. 

The respective velocity distributions of the free-floating 
planets warrant special mention. Planets that are free- 
floating, but still observationally associated with the clus- 
ter, have a very similar median velocity to the whole clus- 
ter. One may naively conclude from this that observing a 
similar velocity distribution for free-floating planetary mass 
objects compared to stars would indicate that they represent 
the very low-mass tail of the IMF, when in actual fact they 
have been liberated from a planetary system. Conversely, 
planets liberated from 5 au orbits which are not associated 
with the cluster appear to have very similar velocities to 
stars, whereas planets originally at 30 au have lower veloci- 
ties. So, in principle one could observe different velocities for 
planets on the periphery of star forming regions and erro- 
neously conclude that their formation scenarios were differ- 
ent. Several studies have already identified potential candi- 
dates for free-floating planetary ma ss objects in young star- 
forming regions, e.g., i n a Orionis (IBihain et al.ll2009h . the 
Orion Nebular Cluster (jWeights et al.ll2009l ). and Ophiuchus 
ijHarvev et al.ll2010l l. Thus far no information about the ve- 
locities of these objects is available but once this is the case 
our results suggest that the analysis of free-floating plan- 
ets' velocities should not be used to draw direct conclusions 
about their formation scenarios. 



5 CONCLUSIONS 

Observations suggest that the initial conditions of star clus- 
ters are cool and substructured. Few authors have exam- 
ined the effects of dynamical evolution in clusters on plan- 
ets through direct A''-body simulations, and none have con- 
sidered initially substructured environments. We have con- 
ducted A'^-body simulations of star clusters in which we place 
a single Jupiter-mass planet at 5 au or 30 au around roughly 
half of the stars in the cluster. The remainder of stars have 
a stellar binary companion, as a high binary fraction is 
observed in most star-forming regions. We have kept the 
level of substructure constant, adopting a fractal distribu- 
tion with fractal dimension D = 2.0. This gives a moder- 
ate amount of substructure, although we note that some 
star form ation regions may have even more primordial sub- 
structure llCartwright fc Whitworthll2004l : ISanchez fc Alfarol 
l2009l : ISchmeiall201ll '). 

We have varied the number of objects in the cluster, 
adopting iVobj = 200, iVobj = 1500 or TVobj = 3000. We 



also vary the initial virial ratio, Q, creating 'cool' {Q = 
0.3), virialised (Q — 0.5) and 'warm' (Q = 0.7) clusters. We 
dynamically evolve each cluster for 10 Myr and examine the 
effects of dynamical evolution on the planetary systems. Our 
conclusions can be summarised as follows: 

(i) In our reference case (Q = 0.3), we find that ~10 
per cent of planets at 30 au are liberated from their host 
stars via interactions in the cluster. Of the planets that are 
'preserved', another ~10 per cent have their eccentricities 
altered by more than 0.1 from an initially circular orbit 
(e = 0). A smaller fraction (typically 5 per cent) have their 
orbital separations altered by more than 3 au. The respec- 
tive numbers for planets originally at 5 au is lower, typically 
by a factor of two. 

(ii) The planets with increased separation and ec- 
centricity would be candidates for future direct imaging 
searches as they spend a large fraction of their orbit at 
large distances from their host star. However, many plan- 
ets (25 per cent) in the e > 0.3, a > 50 au phase space are 
not primordial systems, and have formed via exchange in- 
teractions or capture in the cluster. 

(iii) The fractions of 'liberated' and 'disrupted' planets 
depend on the initial virial ratio of the cluster. If the cluster 
is in cool collapse, it reaches a much denser state than a 
cluster that is in virial equilibrium or expanding. However, 
planets in all clusters are affected by disruption due to low- 
N dynamical decay and regions of localised high density in 
the substructure. 

(iv) A significant fraction of the 'disrupted' planets also 
suffer a significant change in the orbital inclination which 
may lead to further dynamical perturbations with other 
planetary bodies in these systems such as secular planet- 
planet interactions or even planet-planet scattering. The for- 
mer process may then even lead to the creation of close-in 
hot Jupiters with retrograde orbits. 

(v) Planets that are captured during the evolution of 
the cluster are indistinguishable from disrupted planets in 
terms of their eccentricity and semi-major axis, but such 
planets are more likely to be on retrograde orbits, with an 
inclination angle > 90°. 

(vi) The distribution of velocities of liberated planets 
which remain associated with the cluster are similar to the 
stellar velocities, irrespective of the planet's original semi- 
major axis. This suggests that planets with similar velocities 
to stars could be mistaken for a member of the low-mass tail 
of the IMF which formed from core collapse, when in fact 
their formation history was very different. Planets that are 
not associated with the cluster tend to have similar velocities 
to the median non-associated stellar velocity if they formed 
at 5 au, whereas planets originally on 30 au orbits have a 
median velocity that is much lower than the non-associated 
stellar velocity. 

In a future paper we will expand this work to study the 
effects of dynamical evolution in clusters in which planets 
have formed in primordial stellar binary systems. 
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